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1  INTRODUCTION'  AND  EQUATIONS  OI;  MOTION 


Recent  papers  on  the  mechanics  and  mathematical  analysis  of  the  formation  of  adiabatic  shear 
bands  (Molinari  and  Clifton  [19S7]  and  T/avaras  1 1987'.)  have  emphasized  the  late-time,  high- 
temperature  behavior  of  the  equations  for  ihermo/visco/plasticily  in  order  to  judge  whether  or  not  a 
band  will  form  from  a  small  inhomogeneity.  In  fact,  although  the  mathematical  arguments 
presented  in  those  papers  arc  undoubtedly  correct,  they  arc  misleading  in  that  computational 
experience  clearly  shows  that  shear  bands  form  at  intermediate  times  somewhat  after  instability  first 
occurs,  hut  when  the  temperature  in  the  hand  is  still  relatively  low.  It  is  the  purpose  of  this  paper 
to  develop  methods  whereby  the  time  of  intense  localization  (not  the  time  at  onset  of  instability) 
can  be  calculated  or  at  least  estimated  and  to  examine  the  causes  of  band  formation  and  the  effects 
of  the  dominant  parameters. 


Consider  simple  shearing  of  a  slab  of  rigid/plastic  material  of  thickness  x  and  width  2H.  The 
motion  is  assumed  to  be  of  the  form 


x  =  X  +  u(Y.t),  y  =  Y.  z  =  Z,  (1.0) 

where  x,y,z  denote  present  position;  X.Y.Z  denote  reference  position;  and  t  is  time.  This  may 
though!  of  as  an  idealized  version  of  the  Kolsky  bar  experiment  on  a  thin-walled  tube  where  x  is 
the  wail  thickness  and  211  is  the  gage  length.  In  nondimcnsional  form,  the  governing  equations  to 
be  studied  arc 


(xs).y  =  0 

(2.1) 

(XO).,  =  (TS)V,v  +  k(Tl},v),y 

(2.2) 

s  =  Kg(ri)(v.jm 

(2.3) 

X*.,  =  M(K\t))SV,y  , 

(2.4) 

where  -1  <  y  <  +1,  0  <  t  <  and  the  commas  denote  partial  differentiation  with  respect  to  the 
indicated  independent  variable.  The  first  equation  represents  balance  of  momentum  in  the  quasi¬ 
static  approximation  where  s  is  the  shear  stress,  from  finite  element  solutions,  Wright  and  Walter 
[1987],  it  is  known  that  this  approximation  is  very  accurate  up  to  moderate  strain  rates.  The 
second  equation  represents  balance  of  energy  where  i)  is  temperature,  measured  from  some 
convenient  reference  temocraturc;  k  is  thermal  conductivity;  v  -  u.,  is  the  panicle  velocity;  and  v,y 
is  the  velocity  gradient  or  plastic  strain  rate  since  elasticity  is  ignored.  The  first  term  on  the  right 


of  Equation  2.2  is  the  heat  source  due  to  plastic  work  and  the  second  accounts  for  heat  conduction. 
The  third  equation  ;.:n os  the  How  stress  in  temis  of  a  work  hardening  parameter,  tc,  a  thermal 
softening  function,  gift),  which  is  assumed  to  he  monoionicaily  decreasing  with  gt(f)  =  1;  and  a  rate 
hardening  factor  where  m  is  the  rate  hardening  exponent.  The  last  equation  is  an  evolution  law  for 
work  hardening  and  simply  states  that  the  rate  of  hardening  is  proportional  to  the  rate  of  plastic 
work  where  the  factor  of  proportiona'ity  depends  on  the  current  values  of  tc  and  t).  The  inclusion 
of  x  in  Equations  2  1  and  2.2  indicates  that  the  balance  laws  have  been  integrated  through  the 
thickness.  It  is  tacitly  assumed  that  the  field  variables  s.  v,  0,  and  tc  do  not  vary  through  the 
thickness  and  that  r  varies  only  slowly  uiLh  y. 

The  relationships  between  the  nondimcnsional  variables  and  the  physical  variables,  indicated  by 
overbars,  are  as  follows: 

y  =  y/ll  t  =  v,t  (3.1) 

s  =  s/s„  =  pc;Vs„  V  =  v/vll  X'  =  K/Ka  (3.2) 

k  =  k/pcl l"y0  where  s„  =  (bYl))mK  .  (3.3) 

O 

In  Equation  3  yn  is  the  nominal,  applied  strain  rate;  s„  is  a  characteristic  stress  given  by 
Equation  3.3;  p  is  density;  c  is  heat  capacity;  tq,  is  a  characteristic  value  for  the  work  hardening 
parameter,  and  b  is  a  characteristic  time  for  rate  hardening.  Note  that  the  characteristic  stress, 
which  includes  the  rate  factor  (try0)m,  has  been  chosen  so  that  s  =  1  when  tc  =  1,  f)  =  0,  and 
v,y  =  1.  It  is  necessary'  to  define  one  more  parameter,  namely 

a  =  -  g^fO)  =  ust/pc,  where  a  =  -g^(0),  (4.0) 

and  the  subscript  denotes  differentiaron  with  respect  to  the  argument  indicated.  The  form  in  which 
the  equations  have  been  written  is  suitable  for  examining  Lite  effect  of  various  small  defects,  such 
as  variations  in  wall  thickness,  variations  in  strength,  or  variations  in  initial  temperature. 

In  all  problems  to  he  considered  below  the  velocity  is  taken  to  he  constant  on  the  boundaries, 
which  are  assumed  tn  tv  insulated  so  that  the  temperature  gradient  vanishes  there. 


2.  EXACT  RESULTS:  NON-HEAT-CONDUCTING  CASE 


Molinari  and  Cl i Ron  1 1 087 ]  have  given  the  exact  solution  to  the  non-hcat-conducting  case 
when  stress  is  prescribed  on  the  boundary  and  an  approximate  solution  for  the  case  when  velocity 
is  prescribed.  Here,  the  exact  solution  for  the  velocity  boundary  condition  will  be  developed. 

With  k  =  0  in  Equation  2.2  the  thickness  x  factors  out,  and  Equations  2.2  and  2.4  may  be 
combined  to  give  a  differential  equation  connecting  k  and  i3; 

v,t  =  M(tC,l3)0,,  .  (5.0) 

This  implies  that  tc  is  only  a  function  of  0  (and  y  in  inhomogeneous  cases)  k  =  ^(0).  In  turn,  this 
allows  Equation  2.3  to  be  written 


s  =  G  (d)(v.yr.  (6.0) 

where  G(0)  =  K(d)g(0),  and  the  explicit  y-depcndcncc  has  been  suppressed.  The  problem  has  been 
scaled  so  that  G(0)  =  1.  It  will  be  assumed  that  G  has  at  most  one  local  maximum  for  i3  >  0. 
Typically  G  will  increase  initially,  reach  a  maximum,  and  then  decrease  as  shown  in  Figure  lb. 
Combining  Equations  2.2  and  2.3  gives 

=  (xs)fl+m>/m  G'1/m  T'(,+ro)/m  (7.0) 

Because  of  Equation  2.1,  the  term  (Ts)(l+m)/m  is  a  function  only  of  time,  so  that  a  new  time  scale 
T(t)  can  be  defined  by  the  ODE 


dT/dt  =  (is)<1+m>/m  T(0)  =  0.  (8.0) 

Since  the  right  hand  side  of  Equation  8  is  positive,  T  will  increase  monotonically.  However,  when 
a  shear  band  forms,  the  stress  decreases  rapidly,  and  since  (l+m)/m  is  large  when  m  is  small,  T 
will  increase  extremely  slowly  once  the  band  is  well  formed.  Now  the  energy  equation  may  be 
written 


qI/hi  ^  _  ^(1‘mlAn 


(9.0) 


This  has  the  solution 


H(fl)  =  T  t(1+m)/m  +  H(fl0(y)).  (10.1) 

where  H(f»  =  G,/m  dfl.  (10.2) 


Since  G  is  positive,  H  is  monotonically  increasing  and,  thus,  invertible.  Therefore,  i3  is  now 
known  implicitly  as  a  function  of  T  and  y. 


=  §  (T.y).  (11.0) 

From  the  velocity  boundary  condition  and  the  constitutive  relation  in  Equation  2.3,  the  traction  can 
be  found  as  a  function  of  T 

fl  v,y  dy  =  1  =  jlo  (xs)1/m(xG)1/“  dy,  (12.0) 


or  since  xs  is  independent  of  y, 


(xs)‘1An  =  (xCy1/m  dy,  (13.0) 

which  gives  xs  as  a  function  of  T.  Finally,  Equation  13  may  be  used  in  Equation  8,  and  the 
solution  of  a  simple  ODE  by  quadrature  completes  the  problem. 

In  summary,  then,  once  G(O)  has  been  found  from  the  solution  of  Equation  5  and  the 
definition  G  =  &  g,  the  solution  in  parametric  form  is  obtained  by  performing  the  following 
quadratures  in  sequence: 

H(O)  =  H(00(y))  +  T  x(I+m)/m;  (14.1) 


(Ts)'1An  =  (xC)  m  dy; 


(14.2) 


t  =  J*Qr  (xs)-(I+m)An  dT. 


(14.3) 


Note  that  Equation  14  does  reduce  to  the  homogeneous  solution  when  f>0  =  0  and  x  =  1,  since  then 
with  the  aid  of  Equations  14.1  and  14.2  and  the  definition  of  H,  Equation  14.3  reduces  to 
t  =  G  1  dO,  which  can  be  obtained  directly  from  Equation  2.  Also  v,y  =  1  and  s  =  GCd)  in 
this  case. 
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That  Equation  14  gives  localized  solutions  may  easily  be  seen  by  the  following  qualitative 
argument  for  geometric  or  initial  temperature  imperfections.  In  Figure  la  for  a  non-work-hardening 
material,  an  arbitrary  softening  function  g  and  g  raised  to  the  power  1/m  arc  sketched.  Since  m  is 
usually  a  small  number  (0.025  or  smaller  for  many  metals),  the  second  curve  drops  rapidly  toward 
zero.  As  a  consequence,  the  integral  of  g1/m,  also  shown  in  Figure  la,  starts  at  zero,  then  rises 
rapidly,  and  Finally  levels  out  into  a  plateau,  which  has  been  normalized  to  1  in  the  Figure.  The 
situation  for  a  work-hardening  material,  shown  in  Figure  lb,  is  similar,  except  that  in  that  case,  the 
integral  of  G1/m  rises  slowly  at  First,  then  increases  very  rapidly  near  the  maximum  of  G,  and 
Finally  goes  into  a  plateau.  In  the  Figure,  all  curves  have  been  renormalized  as  required  so  that 
each  has  1  as  an  upper  bound. 

First,  suppose  that  z  =  1  and  d0  is  some  nonconstant  function  of  y.  Then  the  initial 
distribution  of  temperature  corresponds  to  some  initial  distribution  of  H,  as  shown  by  the  dark  lines 
near  the  origin  in  Figure  la.  As  time  increases,  T(t)  increases  also,  but  Equation  14.1  shows  that 
for  each  value  of  y,  H  simply  translates  by  the  amount  T,  so  that  at  some  later  time,  H  has  the 
distribution  shown  by  the  dark  lines  near  0.5  and  later  still  by  the  dark  line  near  1.0.  In  turn, 
these  correspond  to  new  distributions  of  -0  as  shown  in  the  Figure.  Since  the  curve  H(i3)  turns  over 
to  the  right,  the  0  distribution  expands,  and  when  the  maximum  value  of  H  encounters  the  plateau, 
the  distribution  of  £  expands  wildly  with  the  maximum  value  of  -0  occurring  at  the  same  location 
as  initially.  In  general,  the  more  sharply  H  turns  onto  the  plateau,  the  more  sharply  will  the 
maximum  in  rise,  and  the  more  peaked  will  the  distribution  in  ■&  become. 

The  graphical  interpretation  for  the  work-hardening  case,  shown  in  Figure  lb,  is  identical. 

Now,  however,  because  H($)  first  turns  up,  the  temperature  distribution  First  becomes  more 
compressed  than  initially.  This  has  been  noted  by  Bums  [1989]  in  a  completely  different  analysis. 
After  G  passes  its  maximum,  and  the  H  curve  begins  to  turn  to  the  right,  the  temperature 
distribution  can  expand,  and  Finally  it  will  peak  sharply,  just  as  before,  when  lire  maximum 
temperature  moves  onto  the  plateau. 

The  graphical  interpretation  is  slightly  different  for  mechanical  imperfection  with  zero  initial 
temperature.  In  this  case,  with  x  nearly  equal  to  1  but  varying  slightly  with  y,  the  H  distribution 
begins  at  zero,  as  indicated  by  the  dots  at  the  origin  in  Figure  lb,  but  according  to  Equation  14.1, 
it  expands  linearly  with  T  as  the  whole  pattern  moves  up  the  H  axis.  Correspondingly,  the  i3 
distribution  starts  at  zero  and  spreads  very  slowly,  as  indicated  in  Figure  lb.  The  process  continues 
until  the  maximum  value  of  H  reaches  the  plateau,  which  occurs  at  the  location  of  the  initial 
minimum  in  d.  As  with  the  case  of  an  initial  temperature  imperfection,  the  temperature  distribution 
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is  compressed  where  H(d)  turns  up  and  expands  and  peaks  up  sharply  when  HO!))  reaches  the 
plateau. 

3.  ASYMPTOTIC  RESULTS:  THE  EFFECT  OF  HEAT  CONDUCTION 

The  result  outlined  above,  which  is  exact  if  the  thermal  conductivity  is  zero,  gives  only  a 
minimum  value  for  the  critical  strain  when  there  is  finite  conductivity.  In  fact,  it  is  more  accurate 
to  refer  to  a  critical  temperature  than  to  a  critical  strain.  Criticality  is  reached  when  the  maximum 
temperature  is  large  enough  to  force  the  H  function  onto  the  plateau.  Since  the  maximum 
temperature  generally  occurs  in  a  narrow  region  with  large  temperature  gradients  on  either  side,  heat 
conduction  has  the  effect  of  delaying  the  critical  condition  by  chopping  off  the  peak  and  spreading 
the  thermal  energy  to  the  sides.  This  effect  was  clearly  seen  in  the  finite  element  calculations 
reported  by  Wright  and  Walter  [1987]  and  plotted  in  their  Figure  8,  rcplotted  here  in  Figure  2. 

The  left-hand  side  of  the  figure  rises  above  the  minimum  as  the  nominal  strain  rate  decreases. 

Since  the  nondimcnsional  thermal  conductivity  varies  inversely  with  the  nominal  strain  rate 
(see  Equation  3),  this  corresponds  to  a  relative  increase  in  the  effect  of  heat  conduction.  The  fact 
that  the  right-hand  side  of  the  figure  rises  above  the  minimum  with  increasing  strain  rate  is  thought 
to  be  an  inertial  effect  and  is  still  not  understood  in  detail.  Since  this  paper  deals  only  with  the 
quasi-static  approximation,  the  right  half  of  the  curve  will  not  be  discussed  further  here,  and  all 
results  obtained  below  will  apply  only  to  the  left-hand  side. 

Only  the  rigid,  perfectly  plastic  ease  will  be  examined,  but  the  defect  that  sets  off  the 
localization  may  occur  in  the  thickness,  the  initial  strength,  or  the  initial  temperature  distribution, 

x  =  1  +  5(y)  ,  k  =  1  +  X(y)  ,  £0  =  i30(y),  (15.0) 

where  max  (I5I,IXI,It30I)  «  1.  In  the  eases  to  be  considered,  there  is  no  work  hardening  so  k  does 
not  evolve.  Consequently,  only  Equations  2.1,  2.2,  and  2.3  apply. 

The  analysis  in  Section  2  introduced  a  second  time  scale,  T(t),  through  the  differential 
transformation  of  Equation  8.  Since  it  turned  out  that  T  was  the  dominant  time  scale  in  the 
absence  of  heat  conduction,  which  is  regarded  as  only  a  small  perturbing  effect  in  the  present  case, 
the  same  transformation  will  be  used  here.  Equation  2  becomes 
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Figure  2.  Strain  at  severe  localization  for  the  initial  temperature  defect  0O  =  0.1(l-y2)9  e'5y2  as  computed  by  finite 
Wright  and  Walter  [1987],  and  by  the  approximate  method  using  Equations  25  and  17. 


(16.1) 

(16.2) 

(16.3) 


(TS),y  =  0 

s  =  Kg(i3)(v,y)m 

xd,x  -  (XKg)'1An  +  k  (TG,y),y, 


where  as  before 

(17.1) 

(17.2) 


dT  =  (xs)<1+m)/m 

w 

(xs)'1/m  =  J’  (xKg)1/m  dy. 


Next  transform  the  dependent  variable  using  Equation  10.2  and  rewrite  the  energy  Equation  16.3  in 
terms  of  H. 


xH,t  =  (xK)-,An  +  k 


dt 

dT 


[(xH,y),y 


H” 

IF 


h.M 


(18.0) 


Upon  ignoring  the  last  term,  which  is  nonlinear  in  H,  the  equation  may  be  written 


or 

x(H-(T),t 

x(H-fT)„ 

=  k  [x(H-fT),y  +  xf,yT],y, 

dT 

=  k  [x(H-fT),y  +  xf,yT],y, 

09.1) 

(19.2) 

where 

(19.3) 

Equation  19.2  has  the  solution 

H  =  ff  +  p. 

(20.0) 

where  P  satisfies  an 

inhomogeneous  heat  equation  with  a  source, 

xp,t  =  k(xp,y),y  +  k(xf,y),yT, 

(21.1) 

P(y,0)  =  H(O0),  p,y  =  -  f,yT  on  the  boundaries.  (21.2) 

Since  X,  5,  and  i}0  in  Equation  15  have  all  been  assumed  to  be  small,  the  initial  value  of  P  is  also 
small,  and  as  a  good  approximation,  Equation  21  may  be  replaced  by 

p„  =  k(p,yy  +  Tf,yy).  (22.0) 
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It  is  convenient  to  split  p  into  two  unknown  functions.  To  that  end,  let  P  =  (4>  -  f)T  +  r,  where  4> 
satisfies 

4>.,  -  1-2.  yy  (23.1) 

4>(y,0)  =  f,  4»,y  =  0  on  the  boundaries,  (23.2) 

and  f  is  def  ined  as  in  Equation  19.3.  Then  the  residual  r  must  satisfy 

r„  =  krlyy +dT(f- «)  (24.1) 

r(y.O)  =  H(f)0(y))  and  r.y  =  0  on  the  boundaries.  (24.2) 


Since  f  -  <J>  is  zero  initially  and  ^  tends  to  zero  as  t  increases,  the  source  term  in  Equation  24  is 
very  weak.  The  Fourier  components  of  r  have  the  representation 


r„  =  f„  c 


k(rtn)' 


'  p  (ck(n’',v-i)  dt’  +  r0  e'1 
J  0  (It  o 


At  early  times,  r  =  0(kt2)  for  all  n  since  ^  is  1  initially.  If  ^  is  bounded  by  e'^  for  some  |d, 
then  r„  ->  0  as  t  — »  Thus,  rn  is  a  small  term.  Finally,  H  has  the  representation 


H(t3)  =  <J)T  +  r  , 


(25.0) 


where  2  and  r  satisfy  Equations  23  and  24.  The  problem  is  completed  by  first  performing  the 
integration  in  Equation  17.2,  which  yields  a  function  of  T  and  t,  and  then  solving  the  ODE  in 
Equation  17.1. 

As  a  first  example,  consider  the  same  problem  considered  by  Wright  and  Walter  [1987]  with 

2 

results  plotted  in  Figure  2.  That  is  let  k  =  1  and  d0(y)  =  0.1(l-y2)9e'5y  .  Then  <j>  =  f  =  1,  and  the 
source  term  in  Equation  24.1  vanishes.  Figure  3  shows  the  calculated  response  for  a  nominal  strain 
rate  of  50s'1.  This  result  is  typical.  Stress  follows  nearly  parallel  to  the  unperturbed  response 
initially,  but  then  falls  off  as  the  temperature  and  strain  rate  in  the  shear  band  accelerate  sharply, 
the  strain  at  intense  localization  is  taken  at  the  point  of  maximum  rate  of  increase  of  strain  rate  in 
the  center  of  the  band.  This  is  the  same  criterion  as  used  previously.  The  solution  that  follows 
from  Equation  25  turns  onto  a  plateau  at  late  times  in  the  same  manner  as  the  full,  finite  element 
solution,  but  the  plateau  occurs  too  soon  so  that  maximum  values  of  temperature  and  strain  rate  are 
underestimated  and  stress  docs  not  drop  far  enough.  The  early  plateau  occurs  because  H”  in 
Equation  18  is  always  negative  and  hence,  the  final  omitted  term  in  Equation  18  is  positive,  which 
would  have  the  effect  of  requiring  larger  negative  curvature  in  the  next  to  last  term  to  balance  it  in 
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the  exact  solution.  Therefore,  the  error  in  the  solution  of  Equation  25  becomes  unacceptably  large 
sometime  before  the  final  plateau  is  reached.  The  criterion  for  critical  strain  stated  above  was  used 
to  obtain  the  approximate  points  plotted  in  Figure  2.  Clearly,  the  approximate  solution  gives  a  very 
good  estimate  of  die  critical  strain  over  several  orders  of  magnitude  of  nominal  strain  rate.  At  the 
lowest  rates,  die  approximate  solution  loses  accuracy  because  heat  conduction  is  not  sufficiently 
accounted  for,  and  at  the  highest  rates,  it  loses  accuracy  because  inertia  terms  arc  ignored. 


As  a  second  example,  consider  a  specimen  with  a  thickness  defect  x  =  1  -  X(l-y2)V5y  . 

Results  for  X  =  0.01  and  X  =  0.001  arc  shown  in  Figure  4.  The  group  of  three  solid  curves  on  the 
left  all  have  the  same  dcf"'t  and  show  that  increasing  nominal  strain  rate  decreases  the  critical 
strain,  which  is  qualitatively  the  same  as  for  temperature  defects.  The  right-hand  curve  shows  that 
a  smaller  defect  increases  the  critical  strain.  Finally,  the  dashed  curve  shows  die  effect  of  ignoring 
the  source  term  in  Equation  24,  which  in  this  case  makes  the  residual  identically  zero.  Figure  5 
shows  the  calculated  residual  for  the  case  X  =  0.01,  %  =  50s'1.  Note  that  the  maximum  is  less  than 
3xl0'\  For  the  same  defect,  but  with  y0  =  500s'1,  the  residual  has  the  same  shape,  but  is  a r  order 
of  magnitude  smaller  everywhere. 


4.  UNIVERSAL  SOLUTION  AND  SCALING  LAW 


Equations  10.2  and  18  provide  the  basis  for  obtaining  considerable  further  insight  into  the 
localization  process  when  m  is  a  small  parameter,  which  is  usually  the  case.  Begin  by  defining  a 
new  variable  z 

l/m  _  I/m 

z  =  -  m  In  f - __  .  (26.0) 

1  -  gl/m 

where  g_  is  the  limit  of  g  as  d  becomes  large.  Then  use  Laplace’s  method  on  Equation  10.2  in 
the  form 


H 


+  0  -  g,/mJ 


Z/rr.]  _(lf> 

J  dz’ 


dz’. 


(27.0) 


Assume  that  dd/dz  can  be  expressed  as  a  power  scries  in  z,  dd/dz  =  a0  +  a:z  +  a2z2  +  ...  where  the 
cocff.Jcnts  a0,  a,,  ...  arc  found  from  Equation  26  by  expressing  g  as  a  function  of  z,  differentiating 
with  respect  to  z  and  evaluating  at  z  =  0. 
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NOMINAL  STRAIN 


O'G  O'Z  O'l  OO  01 

£oix  J  '  ivnaisaa 


*1 


ium  (i-g'  j 

Oo  =  -  =  -  -  =  a 

dz  g’(0) 

d2lX0)  _  1  [  1'^  ,  .An)2  + 

dz2  g’(0)  L  m  m  V  ~  J 


~  a 


(28.1) 


(28.2) 


(28.3) 


When  evaluated  term  by  term,  the  integral  in  Equation  27  gives 

H  =  gi/m  i)  +  (l-gl^XmCI-c^Xao  +  ma,  +  2m2a2  +  ...) 

+  mze'77m(a,  +  2ma2  +  ...)  +  mz2c'2/m  (a2  +  ...) 

+  0  (mz"c'7/m)].  (29.0) 


Since  m  is  small,  the  terms  in  mz"c'z/m  will  be  ignored  at  all  orders  of  n  so  that  a  good 
approximation  for  H  is 


H  -  g,An_  +  01?  (1  -  g1/m),  (30.0) 

where  Equation  26  has  been  used,  and  C  is  a  constant  chosen  so  that  H  -  gi/m  r)  takes  on 
appropriate  values  as  0  tends  toward  infinity  (assuming  that  [g’/m  -  di)  converges).  For 
small  m,  gi^1  may  be  extremely  small.  For  example,  if  g_  =  0.5  and  m  =  0.025,  then  gi^1  =  1012, 
which  for  all  practical  purposes  may  be  set  equal  to  zero.  When  that  is  done,  H  takes  on  the 
simple  form 

H  =  (l-g,ym).  (31.0) 

u. 

With  g  =  exp  (-at)),  this  is  exact,  and  C  =  1.  Plots  of  H,  according  to  Equation  31  or  the  exact 
expression  in  Equation  10.2,  arc  indistinguishable  for  reasonable  choices  for  g  such  as  g  =  1  -  at3 
or  g  =  1/2  (1  +  e  2a\ 

Now,  from  the  definition  of  H  in  Equation  10.2  together  with  Equation  31, 
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H  (0)  =  g,A"  =  1 


and 


H 

mC/a 

ii 


ii"  (0)  =  -  a  (l  -  ) 

tii  C  ^  m(/a/' 


so  dial,  after  a  simple  linear  rescaling. 


A  II  A  T  A  |  A 

li  =  J.  -  ,  T  =  _L-  ,  t  =  __L_  ,  k 
in  (7a  mC/a  mC/a 


mC 


The  basic  ruination  IS  may  be  written  as 


xA.?  =  (tk)  1  m  +  1c 


02.1) 

(32.2) 


(33.0) 


(34.0) 


This  is  a  universal  approximate  equation  for  A  no  matter  what  the  softening  function  g(d)  may  be. 
Since  d?/df  is  given  by  Equation  17.1,  and  from  Equation  31,  it  is  known  that  the  expression 
g  1An  =  1/(1  -  A)  is  a  good  approximation;  Equation  17.2  may  now  be  rewritten  as 


dt  _  mC 

[  r1  (xk)  ,/m 

~a  „ 

dT  a 

1  -  H 

(35.0) 


Equations  34  and  35  imply  that  all  effects  of  the  defect,  whether  in  x,  k,  or  arc  carried  by  the 

A 

universal  function  II. 

Furthermore,  since  Equation  35  shows  that,  for  a  given  defect  in  different  materials,  the  critical 
time  and,  hence,  the  critical  strain,  scales  according  to  the  ratio  mC/a  (with  some  additional  minor 
modification  due  to  rescaled  thermal  conductivity);  it  is  clear  that  the  primary  influence  of  the 
thermal  softening  function  is  through  its  initial  slope  only  (recall  that  a  =  -g'(0)).  The  only  effect 
of  the  exact  shape  of  g(0)  comes  through  the  coefficient  C,  which  may  be  expected  to  have  values 
C  =  1  +  0(m),  according  to  Equation  29. 

Contrary  to  the  emphasis  placed  on  the  behavior  at  large  temperatures  by  Molinari  and  Gifton 
(19871  in  proposing  their  L_  criterion  for  severe  localization,  the  present  analysis  makes  clear  that  it 
is  the  early  behavior  of  the  softening  function  that  determines  when  localization  will  occur.  After 


17 


all,  large  temperature  changes  occur  during  localization,  not  before.  In  fact,  the  temperature 
required  to  carry  H  to  ,99mC/a  (ic.,  99%  of  the  plateau  value)  is  rather  modest  for  reasonable 
softening  functions.  From  Equation  31,  at  this  value  of  H,  g  has  the  value  (.01  )m,  which  decreases 
from  about  .955  to  .891  as  m  increases  from  .01  to  .025.  Then,  for  reasonable  functions  such  as 
1  -  ad,  e1®,  or  1/2  (1  +  c2*0),  all  of  which  have  the  same  initial  slope,  ad  varies  from  about  .045 
to  .047  for  m  =  .01  and  from  about  .109  to  .123  for  m  =  .025.  Notice  that  since  g  is  so  close  to 
1,  m  and  the  initial  slope  of  g  have  more  influence  on  the  corresponding  temperature  than  the  exact 
form  of  the  softening  function.  For  high-strength  steels  (p  =  7800  kg/m3,  c  =  500  J/kg/°K, 
k,,  =  0.5  GPa),  reasonable  values  of  a  are  about  0.1,  which  translates  into  temperature  increases 
when  H  reaches  .99  of  about  45°C  for  m  =  .01  and  about  1 1 5°C  for  m  =  .025.  These  increases  in 
temperature,  reached  just  before  intense  localization,  arc  small  compared  to  those  actually  reached  in 
typical,  fully  formed  shear  bands,  Hartley,  ct  al.,  [1987], 

If  g_  *  0,  the  L_  criterion  of  Molinari  and  Gifton  (1987]  indicates  that  localization  is  not  to 
be  expected  because  H  has  no  limit  as  d  tends  toward  infinity.  However,  Equation  30  shows  that 
this  criterion  is  inadequate  when  m  is  small.  In  fact,  full  finite  element  solutions  for  the  three 
softening  functions  1-ad,  e'4®,  and  1/2  (1  +  e'24®)  all  give  approximately  the  same  time  for  intense 
localization,  Walter  [1989].  The  actual  times  arc  slightly  more  separated  than  predicted  by  the 
scaling  mC/a,  but  the  ordering  is  the  same.  Note  that  the  L„  criterion  predicts  localization  for  the 
first  two  functions,  but  not  for  the  third. 

5.  FURTHER  APPROXIMATIONS  FOR  WEAK  PERTURBATIONS 

Suppose  that  the  wall  thickness  x,  the  strength  k,  and  the  initial  temperature  d0  all  have  a 
weak  cosine  variation  over  the  width  of  the  sheared  section. 

x  =  1  -  5  cosTty 
k  =  1  -  X  costcy 
d0  =  e  cosrcy, 

where  5/m  «  1,  X/m  «  1,  ae/m  «  1.  Then  good  approximations  to  the  solutions  of 
Equations  23  and  24  are 

<>=!+(  Iddn  5  +  x)  cosrcy  (37.1) 


(36.1) 

(36.2) 

(36.3) 
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r  =  1  £  cW  costry,  (37.2) 

A  A 

where  r  has  been  rescaled  in  Equation  25  to  match  H  and  T,  and  the  forcing  function  in 
Equation  24  has  been  ignored  since  it  only  gives  rise  to  a  small  correction.  Now  the  integral  in 
Equation  35  may  be  written  as 

(Tsy1An  =  —1—  f' - - (38.0) 

A  J  0 

1  -  T  (1-a  cos7ty)(l-b  cosJty)(l-c  cost:) 

where  a  =  5/m,  b  =  X/m,  and  c  =  ^ -  I"—  e  +  T  5  +  A.  ^1,  and  a,b,c  <  1  (c<l  can 

L  m  \  m  m/J 

be  justified  a  posteriori).  After  the  change  of  variable  z  =  cosrcy,  Equation  38  is  easily  evaluated 
by  considering  the  related  contour  integral  with  poles  at  z  =  1/a,  1/b,  1/c  and  a  branch  cut  on  -1  < 
z  <  +  1.  The  result  is 


ecs)-1/m 


1  r  a2 

I'  (b-a)(c-a)Vl-az 

+ _ ?! _ ) 

(a-c)(b-c)Vl-c2_  J 


b2 

-  + 

(c-b)(a-b)Vl  b2 


(39.0) 


for  a  *  b  *  c.  To  complete  the  solution  Equation  39  (with  the  above  values  for  a,  b,  and  c 
restored)  replaces  the  integral  in  Equation  35  to  yield  a  nonautonomous  ordinary  differential 

A 

equation  in  t  and  T. 


Rather  than  considering  all  three  types  of  defect  at  once,  it  is  instructive  first  to  consider  them 
one  at  a  time. 


i)  Case  1:  X  =  5  =  0,  e  *  0.  This  is  the  simplest  case.  Equation  35  reduces  to 


4=[o-T>2-(e^c^y](u,n)/2 

dT  L  V  01  '  J 


(40.0) 


In  the  limiting  case  of  k  — >  0  and  m  «  1,  so  that  m  may  be  ignored  in  the  exponent,  Equation  40 
may  be  integrated  exactly  to  yield  an  estimate  for  the  critical  time  (nondimcnsional)  or  critical 
strain 
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=  tCT  = 


mC 


la 


f  l+v'l-(ea/m)2  l 
l  ea/m  J 


mC 

a 


In 


ea/m 


(41.0) 


ii)  Case  2:  8  =  e  =  0,  X  *  0.  Equation  35  reduces  to 
t  (l-c-^) 

dt  ^  \rl  -  (X/mY  (42.0) 

[•^1-t)1  •  (t  4  c'-V  +  ^  -  (Wm)2  tc-]'1  V’" 

-v/o  -  V  •  (t  4  c")’  -* 


As  with  case  1,  in  the  limit  of  1c  — ►  0  and  m  «  1,  so  that  m  may  be  ignored  in  the  final 
exponent,  Equation  42  may  be  evaluated  exactly  to  yield 


Ycr  =  ^  = 


_  Cm 


1 


\T  -  (X/mf 


_  Cm 


In 


2 

77m 


In  1+^1-  (Vm)z  .  Vl+A/m 

VTX/m 


(43.0) 


iii)  Case  3:  X  =  e  =  0,  8  *  0.  This  is  exactly  the  same  as  Case  2  with  X/m  replaced  by 
8(l+m)/m. 


To  treat  all  three  types  of  defect  simultaneously,  consider  the  following  approximation  for  the 
integral  in  Equation  38. 


/’ 
j  0 


dy 

[l-(a+b+c)cos7ty] 


=  [1  -  (a+b+c)2]'1/2. 


Proceeding  as  before  leads  to  the  estimate 

mC 


Ycr  =  f 


CT  T 


In 


2 

(t  A  + 

l+m 

8  + 

V  m 

m 

m) 

(44.0) 


(45.0) 
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Equations  41,  43,  and  45  may  be  regarded  as  approximate  scaling  laws.  Note  that  the 
dominant  effect  comes  from  the  material  properties  mC/a  and  that  the  defect  strength  ae/m,  A/m,  or 
5(l+m)/m  lias  only  a  logarithmic  effect. 

A 

In  evaluating  Equation  40  with  k  =  0,  it  is  clear  that  tcr  occurs  when  the  terms  in  brackets  on 

A  A 

the  right-hand  side  vanish  and  T  =  1  -  ea/m.  When  k  &  0,  Equation  40  must  be  evaluated 

A  A 

numerically.  It  turns  out  t  now  has  infinite  range,  and  T  is  bounded  by  the  curve 
T  =  !  -  e  A  e  ',  so  it  is  not  as  easy  as  before  to  estimate  the  critical  strain,  which  now  must  be 
found  by  actually  computing  the  stress  history  and  judging  where  stress  collapse  occurs.  The  other 
cases  arc  analogous. 

Figure  6  shows  the  stress  response  as  calculated  from  Equations  42  and  44  with  nonzero 
thermal  conductivity.  The  solid  curves  used  Equation  42,  and  the  dotted  curves  used  Equation  44. 
When  the  defect  is  small  enough,  Equation  44  seems  to  give  a  very  good  approximation.  Also 
shown  by  the  dashed  line  in  the  figure  is  the  curve  calculated  from 

Equations  25,  17.2,  and  17.1.  The  nominal  strain  rate  for  these  curves  is  50s1.  Although  the 
solution  for  a  mechanical  defect  has  yet  to  be  calculated  by  the  full,  finite  element  method,  the 
previous  results  for  a  temperature  defect  give  confidence  that  the  dashed  curves  are  accurate  until 
the  plateau  is  reached.  The  figure  shows  that  Equation  42  overestimates  the  critical  time  at  50s'1, 
and  that  the  error  is  smaller  for  smaller  defects.  The  overestimate  occurs  because  of  the 
approximation  (l-a)!An  =  1-a/m  for  a/m  «  1,  which  is  used  in  the  integrand  of  Equation  35  in 
obtaining  Equation  38.  Finally  the  critical  time,  as  calculated  from  Equation  43,  is  shown  by  the 
vertical  interrupted  lines.  Gcarly,  the  simple  formula  gives  a  very  good  estimate  considering  the 
complexity  of  the  full  calculation. 

Equations  41  or  43  or  45  appear  to  give  lower  bounds,  but  that  is  not  quite  the  case. 

Whereas  ignoring  heat  conduction  and  m,  in  tire  final  exponent  1+m,  makes  an  underestimate  of  the 
critical  strain,  the  approximation  for  (l-a)1/m,  as  noted  above,  tends  toward  overestimation  of  the 
critical  strain.  Thus,  the  approximations  tend  to  compensate  each  other  to  some  extent,  although  for 
very  small  defects,  the  last  approximation  can  be  very  accurate. 

6.  CONCLUSION 

The  approximate  methods  used  in  this  paper  reduce  a  complex  calculation  on  a  main  frame  or 
mini-computer  to  a  relatively  simple  calculation  on  a  micro-computer,  since  integration  of  the 
system  in  Equations  17,  23,  24,  and  25  is  easily  accomplished  with  commercially  available  software 
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packages.  The  principal  conclusion  to  be  drawn  is  that  localization  depends  strongly  on  the  low- 
temperature  thermal  softening  and  on  the  rate  sensitivity.  In  fact,  because  of  its  dominant  effect  on 
scaling,  the  ratio  m/a  should  perhaps  be  regarded  as  a  figure  of  merit  that  ranks  perfectly/plastic 
materials  according  to  their  tendency  to  form  adiabatic  shear  bands.  The  shape  characteristic  C  has 
been  omitted  because  its  value  is  always  near  to  one  and  would  be  difficult  to  determine 
experimentally  in  any  case.  The  reciprocal  of  the  ratio  may  be  called  "shear  band  susceptibility." 

Xsb  =  a/m.  (46.0) 

The  susceptibility  controls  not  only  the  critical  strain,  as  has  been  demonstrated  above,  but  in 
an  earlier  paper,  Wright  and  Walter  [1987],  it  was  found  to  play  a  role  in  the  criterion  for  absolute 
stability  and  in  the  early  growth  rate  of  unstable  perturbations.  Susceptibility  does  not  tell  the 
whole  story,  however,  since  the  critical  strain  is  also  dependent  on  many  other  factors;  among  these 
are  defect  size  and  shape  and  thermal  conductivity,  whose  effects  have  been  considered  above,  and 
inertia,  whose  effect  is  yet  to  be  explored. 
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